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Abstract 

In this paper we investigate the dissipativity property of a certain class of Duhem 
hysteresis operator, which has clockwise (CW) input-output (I/O) behavior. In 
particular, we provide sufficient conditions on the Duhem operator such that 
it is CW and propose an explicit construction of the corresponding function 
satisfying dissipation inequality of CW systems. The result is used to analyze 
the stability of a second order system with hysteretic friction which is described 
by a Dahl model. 
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1. Introduction 

Hysteresis is a common nonlinear phenomena that is present in diverse phys- 
ical systems, such as piezo-actuator, ferromagnetic material and mechanical 
systems. From the perspective of input-output behavior, the hysteretic phe- 
nomena can be characterized into counterclockwise (CCW) input-output (I/O) 



dynamics [l|, clockwise (CW) I/O dynamics [21], or even more complex I/O 
map (such as, butterfly map [3]). For example, backlash operator generates 
CCW I/O dynamics; elastic-plastic operator generates CW I/O dynamics and 
Preisach operator can have either CCW or CW I/O dynamics depen ding on the 
weight of the hysterons which are used in the Preisach model [!, [HI, Il6| | . 

In the recent work by Angeli [l|, the counterclockwise (CCW) I/O dynamics 
of a single-input single-output system is characterized by the following inequality 

liminf / y(t)u(t)dt > -oo, (1) 
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where u is the input signal and y is the corresponding output signal. It is 
assumed that u £ U where U is the set of input signals for which y exists and 
is well defined for all positive time. Compare with the classical definition of 
passivity [30j |. it can be interpreted as the system is passive from the input u 
to the time derivative of the corresponding output y. In particular, ((T|) holds if 
there exists a function H : K 2 — > R + such that 

< y(t)u(t). (2) 

Indeed, integrating ([2]) from to oo we obtain (JTJ) . 

Correspondingly, clockwise (CW) I/O dynamics can be described by the 
following dissipation inequality 

liminf / ii(t)y (t) dt > -co. (3) 

The notions of counterclockwise (CCW) I/O and clockwise (CW) I/O are also 
discussed in (20j . 

In our previous results in [13| . we show that for a certain class of Duhem 
hysteresis operator $ : u i— > &(u, yo) := y, we can construct a function Hq : 
M 2 — >■ M + which satisfies 

dH Q {y{t)Mt)) < . (t)tt(t)< (4) 

This inequality immediately implies that such Duhem hysteresis operator is 
dissipative with respect to the supply rate y(t)u(t) and has CCW input-output 
dynamics. The symbol O in Hq indicates the counterclockwise behavior of $. 

In this paper, as a dual extension to we focus on the clockwise (CW) 
hysteresis operator where the supply rate is given by iiy which is dual to the 
supply rate uy considered in [13|. This is motivated by the friction induced 
hysteresis phenomenon in the mechanical system which has CW I/O behavior 
from the input relative displacement to the output friction force. One may 
intuitively consider to reverse the input-output relation of the CW hysteresis 
operator for getting the CCW I/O behavior in the reverse I/O setting. However, 
this consideration has two drawbacks: 1). the reverse input-output pair may not 
be physically realizable (this is related to the causality problem in the port-based 
modeling, such as, the bond graph modeling framework Q); 2). the operator 
itself may not be invertible (for example, if the output of the hysteresis operator 
can be saturated). 

In Theorem [T] we provide sufficient conditions on the underlying functions 
/i and /2 of the Duhem operator, such that it has CW I/O dynamics. Roughly 
speaking, the functions f± and /2 (as defined later in Section 2) determine two 
possible different directions (y, u) depending on whether the input u is increasing 
or decreasing. By evaluating these two functions on two disjoint domains (which 
are separated by an anhysteresis curve), we can determine whether it has CW 
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I/O dynamics using Theorem [T] This is shown by constructing a function Hq : 
R 2 — > R + such that the following inequality 

AH M t)Mt)) < y{t)m (5) 

holds. The function H(j can also be related to the concept of available storage 
function from (30j where, instead of using the standard supply rate yu, we use 
the CW supply rate yu as shown in Proposition 1 in this paper. 

The dissipativity property ([5]) can be further used in the stability analysis 
of the systems with CW hysteresis, such as, a second-order mechanical system 
with hysteretic friction as discussed in Section 4.2. As an illustrative example 
on the application of ([5]), let us consider a mechanical system described by 

VCIX — F -^friction; 
^friction =$(x,y ), 

with the hysteresis operator $ satisfying the Dahl model as follows 

^friction = p ( 1 ~j=T~ ^ max{0, x} + p fl ^ tlon \ min{0, x}, 

where m refers to the mass, x refers to the displacement, F is the applied force, 
p > describes the stiffness constant, Fc > represents the Coulomb friction 
constant and yo is the initial condition of the Dahl model (see, for example, 
(2p||). By taking x\ — x,X2 — x and X3 — Ff r i C ti on as the state variables, we 
can rewrite this hysteretic system into state-space form as follows 

±1 = X 2 , 

"7 m ' \ / \ 

i 3 = P [I - t%) max{0, x 2 } + p [1 + j%) min{0, x 2 }- 

In Section 4.1, we obtain the function satisfying ([S]) explicitly and it is 
parameterized by p and Fc- Using V(x-y, x 2 , x$) = \mx\ + H(j(x^,Xi) as a 
Lyapunov function we have 



V = m±2X2 



dH (x 3l x 1 ) 



dt 

dH (x 3 ,X!) 



= -X3X2 + Fx 2 , 

dt 

< Fx 2 . 

This inequality establishes that the closed loop system is passive from the ap- 
plied force F to the velocity X2- Thus a simple prepositional feedback F — —dx2, 
where d > 0, can guarantee the asymptotic convergence of the velocity X2 to 
zero without having to know precisely the parameters p and Fc- 
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2. Duhem operator and clockwise hysteresis operators 

Denote C 1 (K + ) the space of continuously differentiable functions / : R + — > 
K and AC(R + ) the space of absolutely continuous functions / : R + — > K. Define 

dz(t) ._ z(t+h)-z{t) 

ST — llm h\a+ n • 

The Duhem operator $ : AC(R+) x M — s> AC(R+), (u,y ) h> $(u,y ) =: y 
is described by 0, 0, 0] 

y{t) = fi{y{t),u(t))u+{t) + / 2 (y(t),tt(t))«_(t), y(0) = y o , (6) 

where := max{0,w(i)}, ii-(t) := min{0, The functions /i and / 2 

are assumed to be C 1 . 

The existence of solutions to ((6]) has been reviewed in [18]. In particular, if 
for every £ £ K, /i and / 2 satisfy 

(^i-05)[/i(£ri,0-/i(o5,0] < Ai(0(^-^) 2 , (7) 

(tri-* 2 )[/ 2 (<7i,O-/2(0S,O] ^ - A 2(0(ai -a 2 ) 2 , 

for all o~x, <j 2 € K, where Ai and A 2 are nonnegative, then the solution to ([6]) 
exist and $ maps AC(S. + ) x R — > J 4C(R + ). We will assume throughout the 
paper that the solution to ((6]) exists for all u G AC(M. + ) and yo £ R. 

As a dual definition to counterclockwise (CCW) I/O behavior [![, we define 
the clockwise (CW) I/O dynamics as follows 

Definition 1. An operator Q is clockwise (CW) if for every u G U with the 
corresponding output map y := Qu, where U is the space of input signals such 
that y is well-defined for all positive time, the following inequality holds 

liminf f y(t)u(t)dt > -oo. (8) 

For the Duhem operator inequality (|SJ) holds if there exists a function 
Hcj : M 2 — >• M+ such that for every u € -AC(R+) and j/o G K, the inequality 

dg ° ( ^' M(t)) < y(t)u(t), 0) 

holds for all t where y :— yo). 

In the following subsections, we describe several well-known hysteresis op- 
erators which generate clockwise I/O dynamics and we recast these operators 
into the Duhem operator as in ([5]). 

2.1. Dahl model 

The Dahl model 0, Hi} is commonly used in mechanical systems, which 
represents the friction force with respect to the relative displacement between 
two surfaces in contact. The general representation of the Dahl model is given 

by 



m = p 



^sgn(«(t)) 



^n(l-^sgn(«(t)) )«(/). (10) 
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where y denotes the friction force, u denotes the relative displacement, F c > 
denotes the Coulomb friction force, p > denotes the rest stiffness and r > 1 is 
a parameter that determines the shape of the hysteresis loops. 

The Dahl model can be described by the Duhem hysteresis operator © with 



1 



G 

F~c 

(7 



sgn 



1 



sgn 1 



IT 
(7 



(11) 



(12) 



In Figure [TJ we illustrate the behavior of the Dahl model where F c = 0.75, 
p = 1.5 and r = 3. 




2.2. Bouc-Wen model 

The Bouc-Wen model [25|, [29| is commonly used to model the elastic stress- 
strain relationships in structures. Moreover, it is also used to represent the 
magnetorheological behavior in the MR damper [8] . The general representation 
of the Bouc-Wen model is given by 

y(t) = au(t) - m)WT ~ 7|A(*)l»(*)|y(*)| n-1 , 

where u denotes the displacement, y denotes the elastic strain, n > 1 and /3, £ 
are the parameters determine the shape of the hysteresis curve. 

The Bouc-Wen model can be described by the Duhem hysteresis operator 
© with 

f 1 (a,0 = a-p\o-r-(a\a-\ n - 1 , (13) 

h{a,0 = oc-pW\ n + ^\<- 1 - (14) 
In Figure [2] we illustrate the behavior of the Bouc-Wen model where a = 1 , 
/3 = 1, £ = 1 and n = 3. 
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(a) 



(b) 



Figure 2: The input-output dynamcis of the Bouc-Wen model with a = 1, 
j8 = 1, £ = 1 and n — 3. 



3. Main result 

Before stating our main contribution, we need to introduce three functions in 
the following subsections: an anhysteresis function f an , a traversing function w$ 
and an intersecting function A; these functions will play an important role in the 
characterization of dissipativity and in the construction of the storage function. 
These three functions are defined based on the knowledge of fx and fi . Generally 
speaking, the anhysteresis function f an defines the curve where fx = f^ the 
function w$ describes the trajectory of $ when a monotone increasing u or a 
monotone decreasing u is applied from a given point in the hysteresis phase 
plot, and the intersecting function A defines the intersection of the anhysteresis 
function f an and function u;$ from a given point. The anhysteresis function 
fan and the traver sing function u;$ have the same definitions as given in our 
previous results in [12|. 



3.1. Anhysteresis junction 

In order to define the anhysteresis function, we rewrite fx and fa as follows 



where F, G : R 2 -> K. We assume that the implicit function F(a, £) = can 
be represented by an explicit function a = f a n(0 or £ = g an (a). Such function 
fan (or g an ) is called an anhysteresis function and the corresponding graph 
fan(0)\£ ^ ^} ^ s called an anhysteresis curve. Using f an , it can be checked 
that f\ {f an (£),£) = /2 (/an (0)0 holds. Note also that the functions F and G 
in (TTSl) arc defined by 



h(y(t),u(t)) 
f2(y(t),u(t)) 



F(y(t),u(t)) + G(y(t),u(t)), 
-F(y(t),u(t)) + G(y(t),u(t)) 



} 



(15) 



F = 



h-h 



G = 



fx + h 



2 



2 



3.2. Traversing function uj$, 

For every given point (cr, £) G R 2 in the hysteresis phase plot, let cr, £) : 

[£ , oo) — >• M be the solution a; of 

x(r)-x(0 = J T fi(x(X),X) dX x(0 = a Vrefcoo), 

and let w^O, °"> £) : (— oo,£] — >• R be the solution a; of 

x(t)-x(0= £ f2(x(X),X) dX 1(0 =ff VrG(-oo,e 

Using the above definitions, for every point (cr, G R 2 in the hysteresis phase 
plot, the traversing function cr, £) : R — > R is defined by the concatenation 
of w$,2(-,<7,£) and w$,i(-,ct, £): 

V ; 1 w #i i(r,o-,0 Vt G [e,oo). v ' 

We remark that the function uj$(-,<t, £) defines the (unique) hysteresis curve 
where the curve {(r, oj$(t, er, £)) I T G (— 00, £]} is obtained by applying a mono- 
tone decreasing u to $(-,cr) with u(0) = £, limt-^oo u(t) = —00 and, similarly, 
the curve {(t, cj$(t, cr, £)) I r G [£,00)} is obtained by introducing a monotone 
increasing u to <&(-,cr) with u(0) = £ and limt_ ) . 0O w(i) = 00. 

3.3. Intersecting function A 

The intersecting function A describes the intersection between the anhys- 
teresis curve f an and the curve ui$. The function A : R 2 — > R is an intersecting 
function (corresponding to and f an ) if: i) £J$(A(cr, £), cr, £) = / an (A(cr, £)) 
for all (cr, f) £ R 2 and; ii) A(cr, £) < £ whenever cr > / an (£) and A(cr,£) > £ 
otherwise. This implies that the two functions a;$(-, cr, £) and f an (-) intersect at 
a unique point larger or smaller than £ depending on the sign of cr — f a n(0- I n 
our main result, we also need that eLMeMiHW] exists for every solutions (y, u) 
of©. 

In the following lemma we give sufficient conditions for the existence of such 
intersecting function A. 

Lemma 1. Assume that f\ and fi in (|15| be such that f\, fi are C , Moreover, 
assume that fan is strictly increasing and there exists a positive real constant 
e > such that for all (cr, £) G R 2 the following inequality holds 

whenever a > / on (£) , (17) 

whenever a < f a n(0 ■ (18) 

Then there exists an intersecting function A G C 1 (R 2 , R) such that 
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rr ^TV\ — 


fan s^^^i* 000000 ^ 


>^ \ «(t) 
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Figure 3: Illustration of a Duhem operator with /i(er,£) = e o.5(-i.2<r+£) _|_ q gg 
and / 2 (cr,0 = e ' 5 ( 1,2er_ ^ + 0.83 for all (a,£) € R 2 . The anhysteresis curve 
fan(0 — 0.83£ is shown by the thick solid-line. If the current state be given by 
(y(t), u(t)), the traversing function &<{,(■, y{t), u(t)) is depicted by the dashcd-linc 
and the intersecting point A(y(t),u(t)) is shown by the solid circle. 

(1) A(er, £) < £ whenever a > f a n(£) and A(er, £) > £ otherwise. 

(2) c^(A(<7, 0,^,0 = /a«(A(ff,o). as; 

(3) Moreover, for all u 6 C 1 , y : = j/o), ™e /la-ue i/iaf 4iA(y(t), u{t)) exists. 
The proof of Lemma [T] is given in the |Appcndix A| 

Example 1. In order to illustrate these functions, let us consider the Duhem 
operator $ with fi(a,£) = e o.5(-l.2 CT +£) + q.83 and / 2 (o-,£) = e - 5 ( 1 - 2<T ^ + 0.83 
as shown in Figure [3J It can be checked that the anhysteresis function of the 
operator is f a n(Q — 0.83£ and the functions /i and / 2 satisfy the hypotheses in 
Lemma [TJ With a reference to Figure [3l let the current state of $ be given by 
(y(t), u(t)). In this figure, the traversing function y(t), u(t)) is depicted by 
the dashed-line and the anhysteresis function f an is shown by the thick solid- 
line. The point (y(t),u(t)) is located above the anhysteresis curve, i.e., y(t) > 
fan{u(t)). It can be seen from the figure that the intersecting point A(y(t), u(t)) 
(which is shown by the solid circle) is less than u(t), i.e., A(y(t),u(t)) < u(t). 
This shows that the property (1) in Lemma [T] holds. 



3.4- Duhem operator with clockwise hysteresis 

Based on the three functions w$, f an and A, we define Hq : R 2 — > R + as 
follows 

/■A(o-,£) /-A(c7,e) 
H o (<r,0= fan(r)dr- w 9 (r,tr,t)dT. (20) 

Jo Jf 
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Theorem 1. Consider the Duhem hysteresis operator $ defined in © and (|15l) 

with C 1 functions F,G:M. 2 — > K and with the traversing function lu^ and the 
anhysteresis function f an . Suppose that there exists an intersecting function A 
(e.g. the hypotheses in Lemma\l\hold) . Let the following condition holds for all 
(a, f ) in R 2 

(A) F(o~,£) > whenever a < f an (£), and F(o~,£) < otherwise. 

Then for every u € AC(R+) and for everyyo £ R, thefunctiont — > H&(y(t),u(t)) 
with i?Q as m (|20[) and y := <fr(u,yo), *s rzg/ii differ entiable and satisfies ([5]). 
Moreover, if the anhysteresis function f an satisfies / an (0) = 0, </ien 7?o > 
and the Duhem operator is clockwise (CW). 

The proof of Theorem [T] is given in the |Appendix B| 

Remark 1. In addition to the result in Theorem [7J if f\ and fi satisfy the 
hypotheses given in Theorem{l\ then for every u £ AC(R+) and yo G R, the 
function t — > H^(y(t),u(t)) with as in (|20p is left- differentiate and satisfies 



The proof of this claim follows a similar line as that of Theorem [3 

In order to depict the storage function that is constructed in Theorem[TJ 
we recall again the example of the Duhem operator $ in Example Q] where fx = 
e o.5(-i.2j,+«) +0 83 and j 2 = e 0-5(i.2»-«) +0 83) and it is s hown in Figured Based 

on the functions f an and y(t), u(t)) as shown in Figured] and following the 
construction of the storage function iJ^ as in (|20p. the first component on the 
RHS of (|20|) corresponds to the light grey area in Figure [H Correspondingly, 
the second component on the RHS of ([20]) refers to the dark grey area in Figure 
G3 The summation of these two areas gives the storage function H(j for a given 
state (y(t),u(t)) satisfying JSJ) according to Theorem[T] 

The principle of the construction of in (l2Cfl) can be described in words as 
follows. From a given state (j/(0), u(0)), let us define the trajectory that crosses 
the anhysteresis curve at a given time T by applying either a monotonically 
increasing input signal u(t) — u(0) + t or a monotonically decreasing input 
signal u(t) = u(0) — t. Denote this trajectory by y and the intersecting point 
by (y(T),u(T)). Then the storage function H(j is given by the integral of the 
anhysteresis function from to u(T) minus the integral of y from to T. 

Proposition 1. Consider the Duhem operator $ satisfying the hypotheses in 
Theorem [H Moreover, we assume that the anhysteresis function f an satisfies 
fan(0 = for all £ £ R. Then for every yo,UQ € R, the function as in (|20p 
satisfies 



lim 

hyo- 



H (y(t + h), u(t + h)) - H (y(t),u(t)) 



< y(t)u(t). 



h 




sup 

u£AC(SL + ) 
u(0)=« 
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where y :— &{u, Vp). In other words, H(j defines the available storage function 
fas discussed in f3ul] ) where the supply rate is given by yu ( instead of yu as in 

Bd- 

Proof. As given in the first part of the proof of Theorem [TJ we have 

^-H (y(t),u(t))=u(t)y(t)~ [ iu*(r, y (t), u(t))dr. (21) 

at Ju(t) at 

Integrating (|2ip from t = to T, we obtain 

H (T)-H (Q)= [ y(r)u(r)dr- [ [ ^(r,y(t),u(t)drdt, 

JO JO Ju(t) at 

where u* = A(y(t),u(t)) and we have used the shorthand notation of H^(t) := 
H (y(t),u(t)). _ 

By rearranging the terms in this equation, we arrive at 

y(r)u(r)dr = H o (0) - H (T) - [ [ ^-u> 9 (T,y(t),u(t)drdt. (22) 

JO Ju(t) at 

The supremum of the LHS of (f2"2"]) over all possible u and T defines the available 
storage function where the supply rate is yu. Note that this supply rate is a 



particular class of the general supply rate as studied in [27|, |30| • Since the last 
two terms on the RHS of (|22|) is non-positive, we will show that we can define 
u and T such that these two terms are equal to zero, and thus the supremum of 
the LHS of (f2"2"]) is equal to H^(y(0),u(0)), which is equivalent to H(j(y ,u ), 
i.e., -Hq is the available storage function. 

From a given initial condition (yo, uq), let us introduce an input signal u(t) = 
u (T — t) + tA(y ,u ) for all t E [0,T] and u(t) = A(y ,uo) otherwise. This 
means that we have an input signal u which starts from uq, ends at A(j/oj u o) at 
t = T and remains there for all t > T . Together with the corresponding signal 
y = $(u,yo), we have A(y(t),u(tj) = A(yo,«o) for all t, i.e. the intersecting 
point is always the same. Indeed, this follows from the fact that A(y(t),u(t)) 
remains the same along the trajectories that converge to the intersection point 
(w$(u*,yo,wo),w*) where u* = A(y ,u Q ). 

Following the same arguments as in the proof of Theorem [1] (c.f., the argu- 
ments that lead to Eq. (|B.9[1 ). this input signal ensures that the last term on 
the RHS of (|22|) is equal to zero. Since u(T) = A(y ,uo) for all t > T, we also 
have that H(j(y(t),u(t)) = for all t > T, i.e. the second term on the RHS 
of (I22p is zero using such an input signal. Hence Hq as in (|20|) is an available 
storage function. 

□□□ 

The results given in Theorem Q] can be slightly generalized in order to 
incorporate the case when the Duhem hysteresis operator <J> has saturated 
output. Consider the set D C M 2 which contains all relations of $, i.e., 
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TZ yo U ■= {(y(t),u(t)) G R 2 |y = <f>(u,y ),te R+} C L> holds for all u 6 AC(R+) 
and (yo, w(0)) 6 D. For example, the set D for the Dahl model in Section [2] is 
given by D = (—Fc, Fc) x R. Using D, we can generalize Theorem[T]as follows. 

Proposition 2. Consider the Duhem hysteresis operator $ defined in (j6]) and 
(|15p C 1 functions F, G : D — > R and loitt ffte traversing function w$ and 

i/ie anhysteresis function f an - Assume that the anhysteresis curve is in D and 
there exists an intersecting function A (e.g., the hypotheses in Lemma\]\ hold). 
Assume further that the Assumption (A) holds for all (a, £) in D. Then for 
every u G AC(M.+) and (yo,u(0)) G D, the function t — > H&(y(t),u(t)) with £/q 
as in (|20p and y := $(u,?/o) *s riff/li differentiable and satisfies ([5]). Moreover, if 
the anhysteresis function f an satisfies f an (0) = 0, then > and the Duhem 
operator is clockwise (CW). 

The proof follows the same arguments as that of Theorem [1] 



4. Examples 

4.1. The function for the Dahl model 

Recall the Dahl model as defined in Section [2~T1 and consider the case when 
r = 1. In this case, the Dahl model can be reformulated into the Duhem 
operator as in ([5]) with 



h{<r,0=P\l-y), M<r,$ = p(l+y), (23) 

where p > and F c > 0. It is immediate to check that the conditions as given 
in ([7} are satisfied, which means there exists solution for this Duhem operator 
for all positive time. The anhysteresis function of the Dahl model is f a n(0 = 0. 
Calculating the curve we have 

^(r,^), u(i ))-| _ Fc+(yW + Fc)e -(.-«( t) ) r e (_«,, „(*)]. ™ 

From (f2"5)) and (|2"4"1) . it is immediate to see that the pair (y,u) is well-defined in 
D = (—Fc,Fc) x R. The intersecting function A(y(t),u(t)) is given as follows 

Since f\ and fi in (j2"5)l satisfy the Assumption (A) for all (a, £) G D, the result 
in Proposition 2 holds. 

By denoting u*(t) = A(y(t),u(t)), we can compute explicitly the function 
H(j as follows 
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y(t)+F c 


go 

P 


hi 


-F c 
V(t)-F c 



„ , _ / -Fc(u(t) u*(t)) + f (y(t) + F c )(l - ei^'W-W) tf (t) > 

+ y(*)>o 

»(t)<0. 

Indeed, it can be checked that Hq < ii(t)y(t). 

4-2. Stability analysis of a second-order mechanical system with hysteretic fric- 
tion 

Now, let us consider an example of a mechanical system with the Dahl 
friction model given by mx + dx + kx + $(x) = 0, where m > 0, d > 0, k > 0, 
the hysteresis operator $ is given as in (|2"5|) with p > and F c > 0. As discussed 
before, the functions /i and f 2 satisfy the hypotheses of Proposition 2. 

The state space representation of the system is given as follows 

±x = x 2 , 

±2 =-^1-^2-^, 



±3 -p(l-^Jx 2+ +p^l + ^jx 2 _. 

Using V(xi, X2, Xs) — \kx\ + \mx\ + Hq(xs,Xi), where is as in (l26l) and 
satisfies ([5]), a routine calculation shows that 

V < -x 2 x 3 - dxl + x 3%2 
= —dx\. 

Since the relations of the corresponding Dahl operator (i.e. the set ~R*y 0tU := 
{(y{t),u(t))\y = $(u,y )}) is contained in (—F C ,F C ) xl for all y Q G {-F C ,F C ) 
and u £ AC(M. + ), then it implies that x 3 (which is the output of the Dahl op- 
erator) is bounded and lies in the interval (—Fc,Fc). Additionally, we have V 
which is lower bounded and radially unbounded in the first and second argu- 
ments, i.e. V(xi,X2,Xs) —> 00 as || %\ \\ — > 00. Thus V < —dx\ implies that 
the state trajectory (£1,2:2) is bounded. Moreover, using the boundedness of 
{x\,x>i) and the boundedness of x 3 , an application of the Lasalle's invariance 
principle shows that (xi,X2,x 3 ) converges to the largest invariant set where 
X2 = 0. By analyzing the corresponding state equations, this invariant set is 
given by {(x 1 ,x 2 ,x 3 )\kxi = -x 3 , x 2 = 0}. 



5. Conclusion 

In this paper, we have investigated the clockwise I/O dynamics of a class 
of Duhcm hysteresis operator by obtaining sufficient conditions for the Duhcm 
operators to be CW. The CW property is obtained via the construction of a 
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suitable function satisfying the CW dissipation inequality which can be useful for 
studying stability of systems having CW hysteretic element, such as, mechanical 
systems with hysteretic friction. The sufficient conditions for CW I/O dynamics 
incorporates also the knowledge of anhysteresis function which is commonly 
neglected in the literature of hysteretic systems. For systems identification of 
hysteresis systems, the results provide additional characterization of the Duhem 
operators that can be used to restrict the class of the Duhem operators which 
will be fitted with the measurement data. 
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Appendix A. Proof of Lemma [l] 

Proof. The proof is similar to the proof of [l2|, Lemma 3.1]. Consider the con- 
tinuous function ip : R 3 — >• K defined as yo,Uo) = w$(£, yo, wo) — fan(0- 
Consider also Aq and A\ the two subsets of R 3 defined as, 

^0 = {(£,2/0, Wo) G K 3 , 2/0 > fan(u ) , £ < W } , 
A l = {(£, 2/0, «o) 6 R 3 , 2/0 < /tm(wo) , £ > U } . 

Note that the function f an being strictly increasing by assumption, implies that 
these sets are open sets. Also, the function u;$ satisfies 

/2(w*(£,2/ )Wo),0 V(£,2/o,w ) G ^o , 

/i(^3>(£,2/o,w ),0 V(£,2/o,w ) G ^i ■ 

Consequently, w$(£, j/o, wo) is the solution of ordinary differential equations com- 
puted from C 1 vector field. With [l(J Theorem V.3.1], it implies that w$ is a 
C 1 function in Ao U A\. Moreover, the function f an being C 1 implies that the 
function <p is C 1 in Ao U Ax. With (ITTl) and (fTS)) . the function <p satisfies, 

§£(£,2/o,w ) > e^0 , V(£,y ,uo) € A Q U Ai . 

Consequently, 93 is a strictly increasing function in its first argument in the set 
Aq U Ai. This also implies that, 

2/o, "o) < v(w , 2/o, wo) + e(£ - «o) 

V(£,2/o,w ) G A , 
2/0, wo) > ¥>(wo, 2/0, w ) + e(£ - uo) 

V(C,2/o,w ) G Ai . 

Note that if j/o > /<m(wo), then ip(uo,yo,uo) > and consequently there exists 
a unique real number u* such that tp(u*, j/o, uq) — and (m*,j/o,wo) G Aq. On 



(£,2/o, w ) = 
(£,2/o, w ) = 
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the other hand, if yo < f an {uo), then <p(uo, yo, uq) < and consequently there 
exists a unique real number u* such that <p(u* ,yo,uo) = and (u*,yo,uo) € 
A%. Therefore, by denoting A(yo,uo) — u*, by employing the implicit function 
theorem and using the fact that ip is C 1 , it can be shown that A is C 1 . 

□□□ 



Appendix B. Proof of Theorem [T] 



Proof. The proof of Theorem[T]follows the same line as in our previous work [12j . 
In the first part of the proof we will prove that for all t 6 R + , H?j(y(t),u{i)) 
exists and satisfies (|5|). In the second part we show the non-negativeness of 
H (y(t),u(t)). 

To show that exists, let us denote u* := A(y,u). Using the Leibniz 
derivative rule, we have 



d i . . ... d 

-H (y(t)Mt)) = ^ 



A(y(t),u(t)) 



/a«(V)dT - 



A(y(t),«(t)) 



u(t) 



^{T,y(t),u(t))di 



A(y(t), u(t))f an (A(y(t), u(t))) - A{y(t), u(t))w*(A(y(t), u(t)),y(t), u(t)) 

dw*(r,y(t), «(*)),,_ 

'u(t) 



u(t)u]$(u(t),y(t),u(t)) - 



dt 



u{t)y(t) 



— u$(T,y(t),u(t))dT, 



u(t) 



dt 



(B.l) 



where u*{t) = A(y(t),u(t)) and the last equation is due to uj$>{u{t),y(t), u{t)) = 
y(t) and by the hypothesis given in Lemma [TJ The first term in the RHS of 
(jB.ip exists for all t > since satisfies ©. In order to get ([5]), it remains 
to check whether the last term of (|B.1|) exists, is finite and satisfies 



(B.2) 



«*(*) d 

— «»(T,y(t),u(t))dT >0. 



u(i) 



di 



It suffices to show that, for every r E tt*(i)], the following limit 

lim -[w*(T,y(i + e),u(i + e)) - cjJr,y(t),u(t))} (B.3) 
c\,o+ e 

exist and the limit of (IB.3|) is greater or equal to zero when u*(t) > u(t) and 
the limit is less or equal to zero elsewhere. 

For any e > 0, let us introduce the continuous function u t : K — > R by 



w e (r) = w$(t, y(< + e), u(t + e)). 
More precisely, using (|16j) . w e is the unique solution of 



w e (r) 



y(t- 



u(t+e) 



fi(uj e (s), s)ds Vr > u(t + e) 



(B.4) 



(B.5) 



(t+c) 



f 2 (oj e (s) 7 s)ds vV<it(i + e). 



1G 



Note that u>o(t) = w$(r, y(t), u(t)) as in p^|) for all r € M and 



u € (u(t + e)) = y(t + e) Ve6l + . (B.6) 



In order to show the existence of (|B.3|) and the validity of (|B.2[) . we con- 
sider several cases depending on the sign of u{t) and F(y(t),u{t)). It can be 
checked that the hypothesis (A) on F implies that fi{y{t),u(t)) > f2(y(t),u(t)) 
whenever y(t) < f an (u), and fi(y(t),u(t)) < f 2 (y(t),u(t)) otherwise. 

First, we assume that u(t) > and y(t) > f a n(u(t)). In this case, according 
to Lemma[l] we have u*(t) < u(t). Since u(t) > 0, there exists 7 > such that 
t < u(t) < u(s) for all s in (t, t + 7). It follows from (|B.5|) and assumption (A) 
that for every e G (0, 7): 



dw e (tt(s)) 

— j = / 2 w E K s . u s M s ) 

ds 

> / x (o; e (n(s)),n(s)) u(s) Vs G [*,t + e], 



and the function u> satisfies 
ds 



/i(y(s),u(s)) u(s) V.s G [t,i + e]. 



Since the functions e h-> wo(u(t + e)) and e y(i + e) with e G (0, 7] are two C 1 
functions which are solutions of the same locally Lipschitz ODE and with the 
same initial value. By uniqueness of solution, we get u>o(u(t + e)) = y(t + e). 

This together with the fact that u) e (u(t + e)) = y(t + e) and using the com- 
parison principle (in reverse direction), we get that for every e G [0,7): 

uj e (u(s)) < u>o(u(s)) Vs G [t, t + e]. 

Since the two functions w £ (t) and uiq(t) for r G [u*(t), u(t)] are two solutions of 
the same ODE, it follows that0 w e (r) > u)q(t) and we get that if it exists: 

lim -\ujAt) - uj (t)} < W e[u*(t),u(t)]. (B.7) 

e\0+ e 

Then it is clear that 

lim ~[w £ (t)-wo(t)] > Vt G [u(i),u*(i)]. (B.8) 

e\0+ £ 



In the following, we show the existence of the limit given in (|B.7| by computing 



■"-Otherwise there exist n < T2 such that w £ (ti) = o;o(«(ri)) and cj e (r2) > u)q(u(t2)) 
which contradict the uniqueness of the solution of the locally Lipschitz ODE. 
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a bound on the function e H> i[cj e (r) — wo( r )]- Note that for every e £ [0,7], 



|w e (T)-wo(T)| < \y(t + e)-y(t)\ + 



«(*) 



/ 2 (w £ (a),a)da 



/tt(t+e) 

f 2 (u e (s),s) - f 2 (uo(s),s)ds 

(*) 

< |y(t + e )-y(t)|+ / ( + ) |/ 2 (^( S ),s)|d S 
Jtt(t) 

r u(t) 

I/2 (w 6 (a), s) - / 2 (cj (s), a) I da, 

for all t € [u*(t),u(t)]. By the locally Lipschitz property of / 2 and by the 
boundedness of o> e on [t, for all e € [0,7], it can be shown that there exists 
a, such that a is a bound of / 2 on a compact set. Then 

|« e (r)-wo(r)| < \y(t + e)-y(t)\ 

r u(t) 



L \u) e (s) — w (s)| ds + a\u(t + e) - u(t)\ , 

where L is the Lipschitz constant of / 2 on [w m i n ,cj max ] x [t, u(t)] with 

Wmin = min w c (s), 

(c,«)e[0,7]x[r,tt(t)] 

w ma i = max w c (s) . 

(c,.)€[0 )7 ]x[r,«(*)] CV ' 

With Gronwall's lemma, this implies that for every e e [0,7] 

|w e (r) - wo(r)| < exp((u(t) - r)L) [\y(t + e) - y(t)\ + a\u(t + e) - u(t)\ 
for all r e [u*(t),u(t)]. Hence 

lim -\u e (r) - oj (t)\ < exp((u(t) - r)L) f + a] u(t), 

for all r £ u(t)]. Consequently the limit given in (|B.7[) exists. It implies 

that the inequality (|B.2[) holds when £t(f) > and y(t) > f a n{u(t)). 

For the next case, we assume that u(t) > and y(t) < f an (u(t)). Again, 
according to Lemma[T] we have u*(t) > u(t). Since for every e G (0,7] the two 
functions w £ (t) and u/q(t) satisfy the same ODE foidr e [u(t + e), «*(«)], WO 



2 we have for all r £ [u(t + e),u*(t)] 
doj e (-r) 



dr 



A(w«(t),t) 



du)Q (r) 
dr 



/i(wo(t),t) 
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Vr e [u(t + e),u*(t)}, 



lim -[a> e (r) — w (r)] = 0. 



(B.9) 



We can use similar arguments to prove that (|B.2|) is satisfied when ii(t) < 0. 
Finally, when ii(t) = 0, we simply get 



by continuity of the above bound. 

For the second step, we need to show that is non-negative. Consider the 
case when y{t) > f an (u(t)), we have u*(t) < u(t) and w$(r) > / an (r) for all 
t e [u*(t),u(t)] by Lemma[U Since fani?) belongs to the sector [0, oo) for all 
r G M, we have 



H (y(t),u(t)) = /„n(r)dr+ f an (r) -«»(r, y(t), u(t))dr > 0. 



In case when y(t) < f an (u(t)), we can show the non-negativeness of 7? by using 
similar arguments. 



lim -|w e (r) - w (r)| = 0, 





□□□ 
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